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Xj^ • We propose a scenario in which the dark components of the Universe are manifestations of a single 

Q\ • bulk viscous fluid. Using dynamical system methods, a qualitative study of the homogeneous, isotropic 

background scenario is performed in order to determine the phase space of all possible solutions. The 
specific model which we investigate shares similarities with a generalized Chaplygin gas in the background 
i . but is characterized by non-adiabatic pressure perturbations. This model is tested against supernova type 

o 

!— j ■ la and matter power spectrum data. Different from other unified descriptions of dark matter and dark 

energy, the matter power spectrum is well behaved, i.e., there are no instabilities or oscillations on small 
perturbation scales. The model is competitive in comparison with the currently most popular proposals for 
. the description of the cosmological dark sector. 
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INTRODUCTION 

The crossing of observational data from high redshift supernovae of type la (SNe la) [l[ [2], Q], 



cosmic microwave background (CMB) radiation [5(, matter power spectra [6J|, X-rays from clusters 
of galaxies [3, Q|, and weak gravitational lensing strongly suggests that the present Universe 
is dynamically dominated by a dark sector which is responsible for about 96% of its total energy 
content. It is usually assumed that this dark sector has two different components: (i) dark matter, 
which is supposed to consist of weakly interacting massive particles (WIMPS) with zero effective 
pressure and (ii) dark energy, a mysterious entity which is equipped with a negative pressure. Dark 
matter candidates include axions (a particle present in the multiplet of grand unified theories) and 
neutralinos (light particles present in broken supersymmetric models), but none of these particles 
could be detected until now. The most natural dark energy candidate is a cosmological constant 
which arises as the result of a combination of quantum field theory and general relativity. However, 
its theoretical value is between 60-120 orders of magnitude greater than the observed value for the 
dark energy. An alternative to the cosmological constant is a self-interacting scalar field, known 



as quintesssence. For a brief but enlightening review of these and other proposals, see [10J and 
references therein. 

There exists another route of investigations in which dark matter and dark energy are described 
within a one-component model. According to this idea, dark matter and dark energy are just 
"different faces" of a single, exotic fluid. To the best of our knowledge, the first proposal along 

R n n 

this line was the Chaplygin gas in its original and modified forms [11], [12J, 1 1 31 j - However, this 
unified description of dark energy and dark matter, in spite of many attractive features, seemed 
to suffer from a major drawback: it predicted strong small scale oscillations or instabilities in the 



14|. (On the other 



matter power spectrum, in complete disagreement with the observational data 
hand, the observed matter power spectrum corresponds to the baryonic matter distribution which 
does not exhibit strong oscillations so that this point is still controversial.) The apparently 
unrealistic predictions of the unified Chaplygin gas type models are the result of an adiabatic 
perturbation analysis. 



even avoid this problem 
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ras been suggested that non-adiabatic perturbations may alleviate or 



This paper explores to what extent a viscous fluid can provide a unified description of the dark 
sector of the cosmic medium. The general influence of shear and bulk viscosity on the character 
of cosmological evolution has been studied, e.g. in [l{|, in the context of Bianchi type I models. 
Under the conditions of spatial homogeneity and isotropy, a scalar bulk viscous pressure is the 



2 



has 



23] 



only admissible dissipative phenomenon. The cosmological relevance of bulk viscous media 
been investigated in some detail for an inflationary phase in the early universe (see 2CJ, |2l|, [22 
and references therein). However, as was argued in [lfl I24I ]. an effective bulk viscous pressure 
can also play the role of an agent that drives the present acceleration of the Universe. (Notice 
that the possibility of a viscosity dominated late epoch of the Universe with accelerated expansion 
was already mentioned in 25[). For a homogeneous and isotropic universe, the ACDM model and 



the (generalized) Chaplygin gas models can be reproduced as special cases of this imperfect fluid 
description [lfj. Moreover, in a gas dynamical model the existence of an effective bulk pressure 
can be traced back to a non-standard self-interacting force on the particles of the gas. While these 
investigations were performed for the homogeneous and isotropic background dynamics (a study of 
the background dynamics which is similar to the setup of the present paper was recently performed 
in 26(|), a first perturbation theoretical analysis for a unifying viscous fluid description of the dark 
sector was performed in J27]. 

The bulk viscous pressure p v i sc will be described by Eckart's expression 



.23] Pvisc — ^U-^, 

where the (non-negative) quantity £ is the (generally not constant) bulk viscosity coefficient and 
u^jj, is the fluid expansion scalar which in the homogeneous and isotropic background reduces to 
3H, where H = - is the Hubble parameter and a is the scale factor of the Robertson- Walker 
metric. By this assumption we ignore all the problems inherent in Eckart's approach which have 
been discussed and resolved within the Israel-Stewart theory Q, Q (see also Q Q, 13, Q and 



references therein). We expect that for the applications we have in mind here, the differences are 
of minor importance. 

It is obvious that the bulk viscosity contributes with a negative term to the total pressure and 
hence a dissipative fluid seems to be a potential dark energy candidate. However, a cautionary 
remark is necessary here. In traditional non-equilibrium thermodynamics the viscous pressure rep- 
resents a (small) correction to the (positive) equilibrium pressure. This is true both for the Eckart 
and for the Israel-Stewart theories. Here we shall admit the viscous pressure to be the dominat- 
ing part of the pressure. This is clearly beyond the established range of validity of conventional 
non-equilibrium thermodynamics. As already mentioned, non-standard interactions are required 
to support such type of approach [h! 2^ |. Of course, this reflects the circumstance that dark 
energy is anything but a "standard" fluid. (We mention that viscosity has also been suggested 
to have its origin in string landscape 3l|). To successfully describe the transition to a phase of 
accelerated expansion, preceded by a phase of decelerated expansion in which structures can form, 
it is necessary that the viscous pressure is negligible at high redshifts but becomes dominant later 
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on. 

In ref. [3] a one-component bulk viscous model (BV model) of the cosmic medium was inves- 
tigated in which baryons were not taken into account. As far as the SNe la data are concerned, 
the results of this model were similar to those obtained for a generalized Chaplygin gas (GCG) 
model. But while GCG models predict small scale instabilities and oscillations at the perturbative 
level, the corresponding matter power spectrum of the BV model turned out to be well behaved. 
Of course, the results of a model that does not include baryons cannot be seen as conclusive; after 
all, the observed power spectrum describes the distribution of baryonic matter. In the present 
paper, we perform an advanced analysis which properly includes a separately conserved baryon 
component. Thus we establish a more realistic viscous fluid scenario of the cosmic substratum. 

For the quantitative calculations will use a constant coefficient of bulk viscosity. A qualitative 
analysis with the help of dynamical system methods is applied to visualize the space of more 
general cosmological background scenarios. The corresponding phase space reveals interesting 
new features compared with those generally found in similar models (see, for example, [33] and 
references therein). In the present approach there is an entire singular axis. Solutions of the 
desired type are generated, i.e., solutions for which an initial subluminal expansion is followed by 
superluminal expansion. These solutions are characterized by a bulk viscosity coefficient with a 
power law dependence on the energy density, £ = £o/A where £o =const and v < 1/2. 

We test the results of our model both against SNe la data and the observed matter power 
spectrum. Both the results for type la supernovae and for the matter power spectrum show that 
the BV model is competitive with the ACDM model as well as with quintessence and different 
Chaplygin gas type models. In particular, our minimum x 2 value for the SNe la data is similar to 
the x 2 values of those models. Furthermore, the matter power spectrum represents a good fit to the 
corresponding observational data. We argue that the absence of oscillations and instabilities is a 
consequence of the fact that the pressure perturbations in our model are intrinsically non-adiabatic. 

The paper is organized as follows. In Section[TI]we define our model and specify the homogeneous 
and isotropic background dynamics. Section IIIII performs a qualitative analysis using dynamical 
system methods. In section IIV( the SNe la data are used to restrict the values of the physically 
relevant free parameters. In section [V] the matter power spectrum is determined and compared 
with large scale structure observations. In section IVll we present our conclusions. 
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II. BACKGROUND RELATIONS 



A bulk viscous fluid is characterized by an energy density p and a pressure p which has a 
conventional component pp = ftp and a bulk viscosity component p v i sc = —S,(p) u ^j such that 

P = (3p-t{p)u^ . (1) 

On thermodynamical grounds the bulk viscosity coefficient ^(p) is positive, assuring that the vis- 
cosity pushes the effective pressure towards negative values. In fact, the expression ([T]) is the 



28]. As already mentioned, it follows from 



original proposition for a relativistic dissipative process 

the more general Israel-Stewart theory [29. 130I in the limit of a vanishing relaxation time. We shall 
assume this approximation to be valid throughout the paper. 

Let us consider the cosmic medium to consist of a viscous fluid of the type ([1]), which is supposed 
to characterize the dark sector, and of a pressureless fluid that describes the baryon component. 
Hence, the relevant set of equations is 
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iV - \9^R = Sttg/t^ + T b \ , (2) 

rr ;M = , Tr = (p v + Pv )u^-p v g^ , (3) 

p = p v = f3p v - , (4) 

2^ = , Tr = p m u^ . (5) 

The (super)subscripts v and b indicate the viscous and the (baryonic) matter components, respec- 
tively. Since the matter is pressureless, the total pressure p of the cosmic medium coincides with 
the pressure p v of the viscous fluid. For the homogeneous and isotropic background dynamics we 
shall restrict ourselves to the flat Friedmann-Lemaitre-Robertson- Walker (FRLW) metric, 

ds 2 = dt 2 -a 2 (t)(dx 2 + dy 2 + dz 2 ) , (6) 

n 

favored by the CMB anisotropy spectrum [5]. The dynamic equations then are: 



d\ 2 8vrG 



a , 

2 



(Pb + Pv) , (7) 



2- + \ = -8irGp v (8) 

Pv + 3^(f>v + Pv) = , (9) 

Pb + 3-Pb = , (10) 
a 

Pv = PPv ~ S^oPv ) (11) 



where the bulk viscosity coefficient was assumed to have a power law dependence on the energy 
density p v according to £(p v ) = CoPv with £o = const. The dot denotes differentiation with respect 
to the cosmic time. Of course, not all these equations are independent. Equation (fTU|) is decoupled 
and leads to 

p b = P -f • (12) 

Here and in the following, quantities with a subscript refer to the present epoch and we have 
used ao = 1. From now on we will set (3 = 0, thus assuming the dissipative pressure to be the 
dominating contribution. 

With a change of variables, 

dp v dp v da , . 
Pv = -A = l^A =P - a ' (13) 

where the prime means derivative with respect to the scale factor a, and using relation ([7]). the 
conservation equation for the viscous component becomes 

f/ v + -(pv-S~Zafi)=0 ■ (14) 
a \ a J 

This equation can be recast in an integral form: 

=-2* ai^) , vA ; (15) 
e v {e+l) 3(1 -2v) r 2 ' v ' 

khxa , v = ~ ; (16) 



e + 1 2 



We will be mainly interested in the case v ^ I. 

There is a simple, direct relation between p„ and the scale factor a for v = 0. In this particular 
case, we find, 



Pv = Pb< 





„ k 3 




{ 













(i/ = 0) , (18) 



where -B is an integration constant. It can easily be verified that p v — > a 3 when a — > (pressureless 
matter) and p„ — > cte (cosmological constant) when a — > oo. It is expedient to notice that for = 

nnn 

the total energy density p = p v + Pb coincides with the energy density of a specific GCG [11|, 1121, UM ■ 
Generally, a GCG is characterized by an equation of state (E = const > 0) 

EH 

VGCG = ^ ' ( 19 ) 

Pgcg 

6 



which corresponds to an energy density 

F 



PGCG 



E + 



(20) 



a 3(l+a)_ 

k 1/2 , p _ o 1/2 • 



where F is another (non-negative) constant. With the identifications E = and F = -B/O fe g & 

becomes obvious that for a = —1/2 the energy densities p = p v + Pb from ([18]) and pacG m ([20]) 
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coincide. The analogy between dissipative and GCG models was also pointed out in 

The structure of the total energy density p = p v + Pb allows us to perform a decomposition of 
the viscous-baryon system into three non-interacting components plus the baryonic fluid: 

P = Pi + P2 + P3 + Pb , (21) 



with 



rp2 EF F 2 ~ Pbo Pbo fc>0 , 

pi = E , p 2 = 2— , p 3 = - , p b = — , (22 

a A a 6 



a 2 



and 



E = \^fm , F = ^B . (23) 
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The individual equations of state of the different components are: 

Pi = ~Pi , P2 = ~^P2 , P3 = , Pb = . (24) 

Consequently, our system is equivalent to a mixture of two dark energy type fluids, one dark 
matter fluid and the baryon component. While the BV model is a unified description of dark 
matter /energy, the option of a decomposition may be useful in comparing the scenario with some 
observational data. The X-ray measurement of galactic clusters [34] is an example which requires 
a separation of a matter component of the cosmic medium. For a Chaplygin gas, different decom 
positions have been proposed in the literature 
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361 ] . Moreover, the decomposition (|2ip with 



- (|24|) reveals, that the background dynamics of our model differs from that of the ACDM model 
(including a baryon component) by the existence of the second dark energy component with energy 
density p2- 

The relevant cosmological quantities are the present density parameters £I v q and £lbo which 
represent the fractions of the viscous fluid and of the pressureless matter, respectively, with respect 
to the total density (remember that we restrict ourselves to the flat case), the Hubble parameter 
today Hq, the present value of the deceleration parameter qo, and the age of the universe to- The 
apparently new parameter of the present model, the viscosity coefficient, can be expressed in terms 
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of the other quantities. In fact, the relation (|18p can be rewritten in a more convenient way. Using 
Friedmann's equation today, we have 



m 



8ttG 



(PbO + Pvo) 



_ 8ttG 
6H Q 



Q V + ^bO = 1 
8?rG 



n 



v0 



3H$ 



PvO 



Combining (|18p with (|25p and (|26p we obtain the relation 

B + l 1 



60 



(25) 
(26) 

(27) 



Now, if we use equation ([8]) with p = p v = — 3£o§ and the definition of the deceleration parameter 



9o 



2g + 1 



24vrG^o 
In 



From the definition of k in equation (117p we have 

24ttG£o 1 



Using equation ([8j), evaluated today, it follows that 

l-2q 



k 



60 



which leads to 



B 



21 + go 



3 vHw) 

Hence, the viscous fluid energy density becomes, 



Pv = Pb{ f 2 (a) - 1 



1 



2(1 + g ) + (1 - 2g )a : 



3/2 



(28) 



(29) 



(30) 



(31) 



(32) 



3\/Q&q 

All these relations will be useful for defining the observables that will be constrained later on by 
the SNe la and matter power spectrum data. 



III. A DYNAMICAL SYSTEM ANALYSIS 

The relations of the previous section that will be used to compare the theoretical predictions of 
our model with observational data, are valid under the restriction v = 0. It is important to have 
an idea of how serious this restriction is. If the case v = is very particular, our results will not 
be conclusive. Therefore it is desirable to have at least a qualitative analysis of the general case 
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FIG. 1: Phase diagrams for v < 1/2 (upper pane) and for 1/2 < v < 3/2 (lower panel). For v > 3/2 the directions 
of the arrows in the lower panel have to be reversed. 

([7|)- (fTT]) which, because of its complexity, cannot be solved analytically. Such an analysis could 
reveal whether or not the case v = is, in a sense, typical. To this purpose we shall use dynamical 
system techniques for a system of non- linear differential equations 32]. 



The system of equation (jT])- (|lip can be recast in a more convenient form by using geometric 
unities 8ttG = 1 and fixing the scale so that £o = 1. Defining x = ^ and y = p v , the following 
system of equations is obtained: 



--x(x - y v ) = P(x,y) 

-3xy (I - Zxy"- 1 ) = Q{x,y) . 



(33) 
(34) 



Since the density pb = 3x 2 — y of the baryonic matter must be non-negative, the solutions that 
have physical meaning are those which satisfy y < 3x 2 . 

In a first step we have to identify the critical points (xq, yo) for which P(xo,yo) = Q(xo,yo) = 0. 
These points are either attractors or repellers or saddle points. The nature of a critical point is 
determined by the eigenvalues of the matrix 



/ Px(x ,y ) Py{x ,y ) \ 

\Qx(xo,yo) Q y (x ,yo) J 
where P x i = Jj£, x l = x,y and Q x i = J^-. In general, there are two eigenvalues for the 
characteristic equation of this matrix. If both eigenvalues are positive, than the critical point is 
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repeller; if both are negative, the critical point is an attractor; if one is positive and the other one 

n 

is negative, than the critical point is a saddle. With the help of a suitable transformation [32], the 
critical points at infinity can be studied. A projection onto the Poincare sphere [32J will then allow 
us to represent the entire phase space in a finite region in the x — y plane. 

It turns out that the entire Oy axis is a stationary solution for the system, representing a 
static space-time. A singular point, denoted by X, is (xo,yo) = ^S 1 - 2 " , 3 1 - 2 " J where v / i. The 
characteristic equation for this critical point is given by A 2 + ao (| — v) A + &o (2 ~~ u ) = ^ wnere 
ao = 3 1 - 2 " . The eigenvalues are Ai = [y — 2) ®o and A2 = —do- The phase diagram depends on 
the value of the parameter v. If v < 1/2, the critical point X is an attractor: all solutions converge 
to it. This point represents the de Sitter phase. Hence, all solutions approach a de Sitter phase 
asymptotically. 

If 1/2 < v < 3/2, the point X becomes a saddle point and all solutions approach a Minkowski 
space-time in the origin. If v > 3/2, X is still a saddle point, but the directions of the curves change: 
all solutions leave the Minkowski space-time. These different behaviors are shown in figure (pQ). 
For the point X we have 3xq = y$, which corresponds to a vanishing baryon density. (Recall that 
the range y > 3x 2 is unphysical.) The deceleration parameter q can be written as q = — 1 — j^j, 
equivalent to q = — 1 — Combination with Eq. ()33[) shows that there is accelerated expansion 
q < for x < 3y u . The special case v = has (xq, yo) = (1, 3) and accelerated expansion for x < 3. 
This solution is of the type of the curves that approach the point X from right below in the left 
panel of Fig. [TJ Notice that there is a projection effect due to the circumstance that infinity is 
mapped onto one single point on the x axis. The curve that points from the origin towards X in 
the left part of Fig. [T] corresponds to a phantom scenario with pb = which separates the physical 
and unphysical regions. We are not interested here in this type of models. 



IV. CONSTRAINING THE MODEL USING SNE IA DATA 

The type la supernovae allow us to obtain informations about the universe at high redshifts. 
Today, observers have detected about 300 of these objects with redshifts up to almost 2. In what 
follows, we will use the more restricted sample of 182 SNe la of the GoldOG data set which consists 
of objects for which the observational data are of very high quality [37]. The relevant quantity for 
our analysis is the moduli distance fi, which is obtained from the luminosity distance Dl by the 
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-0.493 



TABLE I: The best-fitting parameters, i.e., when \ 2 is minimum, for the BV model, the GCG model, the 
traditional Chaplygin gas model (a = 1) and the ACDM model (a = 0) for a flat Universe (fifco = 0). Hq 
is given in km/Mpc.s, E = -prrp m units of c 2 (c - speed of light) and to in Gy. 



relation 



"= 5ios O +25 



(35) 

The computation of the luminosity distance follows the standard procedure [l[ [3]. First, we 

recall its definition 0,0, 



D L = - = (l + z)r , 



(36) 



where r is the comoving coordinate of the source and z is the redshift, z = — 1 + -. The coordinate 
r can be obtained by considering the propagation of light ds 2 = 0. This implies, 

" a da 



ds 2 = c 2 dt 2 - a 2 dr 2 = 



a a 



From Friedmann's equation one obtains, 



(37) 



(38) 
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TABLE II: The estimated parameters for the BV model, the GCG model, the traditional Chaplygin gas 
model and the ACDM model for a flat Universe (f^o = 0). We use the Bayesian analysis to obtain the peak 
of the one-dimensional marginal probability and the 2 a confidence region for each parameter. H$ is given 
in km/s/Mpc, E = E ^_ p in units of c 2 (c - speed of light) and to in Gy. 

In terms of the redshift z, the luminosity distance may be expressed as 

c f z dz' 

Dl = 3(1 + z) _y o ___ 7 __ , ( 39) 

or, using B and f2 m n, 

c f z dz' 

DL = {l+Z) lfJ, 1 + ^5(^ + 1)3/2-1) ' (40) 

In order to compare the theoretical results with the observational data, the first step is to 
compute the quality of the fitting through the least square fitting quantity x 2 . We adopt the HST 
(Hubble Space Telescope) prior [41J for H$, as well as the cosmic nucleosynthesis prior for the 
baryonic density parameter Qboh 2 (where h is Hq divided by 100 km/s/Mpc), which leads to 

2 _ y Qjjj ~ ^0 (fr-0-72) 2 (Q b0 h 2 - 0.0214) 2 
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The quantities \i§ j are the distance moduli, observationally measured for each supernova of the 
182 Gold06 SNe la dataset 3 

and the /Xq i are the corresponding theoretical values. The a 2 i 
represent the measurement errors and include the dispersion in the distance moduli, resulting from 
the galaxy redshift dispersion which is due to the peculiar velocities of the objects (cf. Q]). 

The smaller the x 2 > the better the agreement between the theoretical model and the observa- 
tional data. Table U shows that the x 2 value for BV model is competitive with the GCG model (cf. 
Eq. ([20]) ). the traditional Chaplygin gas model (cf. Eq. (f20|) with a = 1) and the ACDM model. 
Compared with the ACDM and the (generalized) Chaplygin gas models, the age of the Universe is 
greater for the bulk viscous model and its deceleration parameter is less negative, i.e., the Universe 
is less accelerated. 

Using Bayesian statistics, a probability distribution can be constructed from the \ 2 parameter 



38J. In the present case, there are three free parameters: the Hubble parameter today Hq, the 



pressureless baryonic density parameter Q^o an d the auxiliary quantity B which is connected to 
the present value of the deceleration parameter (cf. (|3ip ). Marginalization over one or two of 
these parameters will lead to corresponding two- or one-dimensional representations. The details 
of the Bayesian statistics and computational analysis, rel ying on BETOCS (BayEsian Tools for 
Observational Cosmology using SNe la), are given in ref. 



The main feature of the parameter estimation for the viscous fluid model is the Gaussian-like 
distribution around the preferred value. 

This can easily be seen both in the table |TT] (almost symmetric error bars) and in the two and 
one-dimensional diagrams of figures [2] and El In figure [2] we show the two-dimensional probability 
distribution for the parameters B, fi&o and Hq at one, two and three sigma, corresponding, re- 
spectively, to 68%, 95% and 99% confidence levels. The corresponding one-dimensional probability 
distributions are shown in figure The main features of the model are: the age of the universe 
is considerably larger than for the ACDM and GCG models, around to ~ 16 Gy, with a small 
dispersion; the baryon density is around 0.05, a little larger than predicted by nucleosynthesis, but 
there is a marginal agreement; the Hubble constant is around 62 km/s/Mpc, i.e. much smaller 
than the one predicted by the CMB, which is around 72 km/s/Mpc, but in agreement with the 
ACDM and GCG models when only the data from SNe la are used; the deceleration parameter is 
bigger (smaller absolute value) than in those models. For the parameter estimations in the ACDM 
and in the GCG models, see ref. [38| and references therein. 

The Gaussian nature of the distribution is responsible for the fact that the best fitting set of 
values is very close to the marginalized parameter estimations, with a small dispersion. Generally, 
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FIG. 2: The two-dimensional plots of the Probability Density Function (PDF) as function for the parameter B, the 
baryon density and the Hubble constant. The joint PDF peak is shown by the large dot, the confidence regions of 1 a 
(68, 27%) by the red dotted line, the 2 a (95, 45%) in blue dashed line and the 3 a (99, 73%) in green dashed-dotted 
line. 

this is neither true for the ACDM model nor for the GCG model. Again: the most remarkable 
differences between the viscous model and the ACDM model are the larger age of the universe and 
the smaller absolute value of the deceleration parameter. 



V. THE MATTER POWER SPECTRUM 



A. The perturbed equations 



In order to write the first order perturbed equations, we re-express the field equations as 



(42) 



JT ;M = , Tr = (p v + Pv )u^-p v g^ , (43) 

Pv = , (44) 

r r ;M = - T£ v = p b u»u v . (45) 

We introduce the quantities, 

9»v =9Ku + h lu/ , p = p + 5p , u>* = u fl + 5u^ , (46) 
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FIG. 3: The one-dimensional plots of the PDF as function for the parameter _B, the baryonic density, the Hubble 
constant and the deceleration parameter. The joint PDF peak is shown by the large dot, the confidence regions of 1 a 
(68, 27%) by the red dotted line, the 2 a (95, 45%) in blue dashed line and the 3 a (99, 73%) in green dashed-dotted 
line. 

where g^v, p and are the known background solutions for the metric, the energy density and 
the four velocity, respectively, and huvi $P an d are the corresponding perturbations. We will 
perform the calculations using the synchronous coordinate condition, 



V) = 



(47) 



Since we will be interested mainly in perturbations that are inside the horizon, the choice of the 
gauge is not really essential. 

The perturbed Ricci tensor takes the form, 



6R 



■flV 



■p Xfip v 



(48) 



where 



Afj,u 



pa 



l <rp\v ~ '*(jv\\i ' L p,v\(j 



(49) 
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are the perturbations of the Christoffel symbols. The relevant non-vanishing components are: 



4, = -\(%) 



X% = ■ (51) 



At j 



~ 2a? { kik,j + hij ' k ~ hij,k } ' ^ 
For the components of the perturbed energy-momentum tensor we have 

5T 00 = 5p , (53) 

5T m = (p + p)5u i , (54) 

ST ij = h ij p-g ij 5p . (55) 



Using the definitions, 



hkk . . ,, • Sp b 8p v 

a 2 Pb £o 



and performing a plane wave decomposition of all perturbed functions according to 

5f(x,t) = 5f(t)e rk - £ , (57) 

where k is the wave vector, we find, after a long but straightforward calculation, the following set 
of first order perturbed equations: 

h + 2% = 8ttG(5p v + 3Sp v ) + 8irGp b 5 b , (58) 

8p v + 3^(Sp v + 6p v ) + (p v + Pv )(&-^j=0 , (59) 
a k 2 

(pv +Pv)@ + (p v +Pv)@ + 5-(p„ +Pv)@ ^5p v = , (60) 



a a* 



h = \ ■ (61) 

Replacing h by using the last relation, we obtain 

S b + 2-4 - 4irGp b 5 b - 4nG8p v - 12TrG£ Ap v = , (62) 
a 

5p v + 3-(Sp v + CoAp„) - (p v + p v )Ap v = , (63) 
a 

a k 2 

(p v +p v )@ + {p v +p v )Q + 5-(p v +p v )Q-—a 2 Ap v = , (64) 

@-8 b + Ap v = . (65) 
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With the help of the constraint (|65|) , this system can be further reduced to a set of three coupled 
differential equations: 



5 b + 2-5 b - ^KGp h 5 h - 4irG5p v - 12irG£ Ap v = 
a 

5p v + 3^(5p v + £oAp v ) - (Pv + Pv)Ap v = 

(Pv +Pv)$b + (Pv +Pv)'k + 5-(pv +Pv)k = 

a 

a k 2 
{pv +Pv)Ap v + (p v + p v )Ap v + 5-(p v +p v )Ap v + —a 2 Ap v 

a so 



(66) 
(67) 

(68) 



Now it is convenient to introduce the background equations and the quantities in Eq. (|32j) from 
which we find 



H 2 
AirGpv 

12irG£ 



3 &bo 
2^3" ' 

2 a a 

5(1 -2go) • 



The following relations will be useful later on: 

1 f(a) 1 



Pv 
Pv 

pv 
Pv 



9{a) 



,3/2 



3 f(a) 



1 y/U 



60 



a 



-3-{p v +p 
a 

1-2^ 



(Pv + Pb+Pv) 



(69) 
(70) 
(71) 

(72) 

(73) 
(74) 



In a next step we divide the set of equations (|66|) - (|68|) by -ffg and redefine H^Ap,, — > Ap„. 
Again, after a fairly long calculation, the perturbed equations reduce to the following system of 
coupled linear differential equations: 



4 + 2 a -5 b - - 3 -^(f 2 - 1)5 V -1(1- 2q )A Pv 

a ' 



2a 3 2a 3 







a 



ii + //i^|2-(i+j) + L r 

a 2 



S„ - 3-5<J„ - (1 + 2g)A Pv = 
a 



( l +//)Ap r + <{ 2^(l + 5 ) + ^° 



+ [kl H 



■ (l-2g )a 
9^o(/ 2 - 1) 



+ 9 



Ap v 



(75) 
(76) 



(77) 
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where Ijj = c/Hq. 

Now, we re-express the perturbed equations in terms of the variable a according to (|13|) . which 
implies 



6 = 6*0, , 5 = a 2 5" + a6' 



(78) 



We also need the following relations, which result from the background equations: 



■ 2 r> f 

a — "60 — 



If / 2 "1 



a a 2 2a\ * p 



+ 1 



With the redefinition A„ = ^jpS the final form for the set of perturbed equations is: 

/ 3 f-l , 3 \ 3 4 

d & + < ! - ofl— ^ — >0 h - - 



2" / 2 a ' 2a J "° 2 / 2 a 2 
3 / 2 -l , 3g / 2 -l 
2 / 2 a 2 2a f " ' 

6' v - 3^6 V - (1 + 2g)A v = ; 

(i+sKM 2«)^+ = 



2a^ ' p 

2 / 



2a^ ' p 
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3pn b0 



A, 



(79) 



(80) 
(81) 

(82) 



B. Numerical integration 



Evidently, equations (|80p - (]82p are too complicated to admit analytical solutions. Hence, we 
proceed integrating these equations numerically. In order to do so, an important problem is to 
fix the initial conditions. A scale invariant primordial spectrum (as predicted by the inflationary 
scenario) is assumed. But, since we are interested in the power spectrum today, we have to follow 
the evolution of this spectrum until the present phase. 

The correspondin g pr ocedure for the ACDM model has been carried out in terms of the BBKS 



transfer function 



43, 



q = q(k) 



441 ] . This function is characterized by 
k 



n° M he- n B-( n B/ Q M) 



(83) 



(nr)Mpc" 1 ' 

where fi^L and O'g are the density parameters of dark matter and baryonic matter of the ACDM 
model. The baryonic transfer function is approximated by by the numerical fit 

In (1 + 2.34q) 



42J: 



T(k) 



2.34 q 



1 + 3.89g + (16.1<?) 2 + (5.46g) 3 + (6.71q) 



z\3 



~\4 



-1/4 



(84) 
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The quantity of interest is the power spectrum for the baryonic matter today. This spectrum is 
given by 

P(k) = \8 b (k)\ 2 = AkT\k)4^\, (85) 

where Qj, is the total density parameter of the ACDM model. The function g in (|85p is the growth 
function 

-l 



0V 7 -0 A+ ( 1 + §)(1 + ^ 



(86) 



43, 



Ml Qrms Of 



in which 17 a represents the fraction of the total energy, contributed by the cosmological constant. 
The normalization coefficient A in (|85p can be fixed by using the COBE measurements of the CMB 
anisotropy spectrum. This coefficient is connected to the quadrupole momentum 
this spectrum by the relation 

A = {2l H f ^ , (87) 

where To = 2.725 ± 0.001 is the present CMB temperature. The quadrupole anisotropy is taken as 

Q rms = l8fiK. (88) 

This value is obtained from the COBE normalization and is consistent with the more recent results 
of the WMAP measurements which use the prior of a scale invariant spectrum [44]. [4^. Taking all 
these estimates into account, one may fix 

A = 6.8 x 10 5 h' 4 Mpc 4 , (89) 

which corresponds to the initial vacuum state for the density perturbations used in the BBKS 
transfer function. 

We recall that all the relations (|83 p -(|89 p are valid for the ACDM model. Strictly speaking, 
we would have to redo all the calculations to obtain the analogue of the expression ([85 p for our 
bulk viscous model. However, since we are interested in the shape of the spectrum, there exists 
a simpler way which relies on the circumstance that the growth function g in (|86|) is a pure 
background quantity. Use of the ACDM growth function and the ACDM initial values also for our 
model will result in an overall shift of the spectrum P without affecting the structure of the latter. 
The resulting data will then be a fit to a "wrong", i.e. unphysical ACDM model, in which the 
composition of the dark sector differs from the standard one with roughly 30% of dark matter and 
70% of a cosmological constant. However, this procedure provides us with the correct structure 
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of the spectrum for our bulk viscous model. Keep in mind that the background dynamics of our 
model differs from that of the ACDM model by the existence of a second dark energy component 
(component 2) in Eq. f)22[) . 

Well inside the matter dominated phase (say, z ~ 500), we use the ACDM initial conditions, 
insert them in equations (|80|) -(|82|). and let this system evolve. Then we compute the power spec- 
trum P(k) = 5?(k) at z = 0, equivalent to a = 1. Finally, we compare the theoretical results 
with the power spectrum obtained by the 2dFGRS observational program. The spectrum depends 
essentially on two parameters: £lbo and go- From now on, relying on the results of the type la SNe 
analysis and the primordial nucleosynthesis constraint, we fix f2&o = 0.04 (according to the type la 
SNe analysis of the previous section, this is marginally satisfied at least at 2 — a). 

In figures U] and [5] we present the best fitting for qo = —0.3,-0.4,-0.5,-0.6. To have a good 
agreement, the ratio between the quantity of dark matter/dark energy in the "wrong" ACDM 
reference model has to be in the range from 1 for qo = —0.3, to about 2/3 for qQ = —0.6. It is seen 
that the baryonic matter power spectrum of the viscous fluid model is in good agreement with the 
observational data for values of qo that are close to those predicted by the type la SNe analysis, 
that is, qo ~ —0.4. We remark en passant that a variation of the value of f^o within the interval 
predicted by the type la SNe at 2 — a does not change substantially the results presented in figures 
H]and[I3 In any case, there is no instability in the computed power spectrum. This is true not only 
for the baryonic fluid power spectrum but also for dark viscous fluid power spectrum. We note, 
however, that as qo approaches zero, there is a power depression at large scales (small k) compared 
with the ACDM reference model. 

C. A physical interpretation 

In the homogeneous and isotropic background, the BV model is equivalent to a GCG with 
a = —1/2. The qualitative differences between both models at the perturbative level can be traced 
back to a difference in the pressures that characterize the cosmic medium in each of these cases. 
In the background both pressures coincide. But while the perturbations for the GCG are entirely 
adiabatic, this is not the case for our present model. 

In the background we have p = p v = —3H^o for the viscous fluid and p = -Ep l l 2 for the GCG 
with a = —1/2 (cf. (|19p ). Use of Friedmann's equation 3H 2 = 8irGp (for the spatially flat case) 
reveals that in both cases we have an equation of state p cx —p 1 / 2 , where the constants are related 
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FIG. 4: Predicted power spectrum for the dark viscous fluid compared with the 2dFGRS observational data for 
go = —0.3 (left) and go = —0.4 (right). The initial conditions using the BBKS transfer function correspond to 
"wrong" ratios (see text) of dark matter and dark energy: QdmO = 0.49 and SIao = 0.51 (left) and QdmO = 0.52 and 
Ovo = 0.48 (right). 
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FIG. 5: Predicted power spectrum for the dark viscous fluid compared with the 2dFGRS observational data for 
go = —0.5 (left) and go = —0.6 (right). The initial conditions using the BBKS transfer function correspond to 
"wrong" ratios of dark matter and dark energy: fldmo = 0.44 and Oao = 0.56 (left) and fldmo = 0.52 and JIao = 0.48 
(right). 
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by E = ^ V24vrG. The first order pressure perturbations for the GCG are 



Sp 




(GCG model) 



(90) 



where 



P 



P 
P 



(GCG model) 



(91) 



P 



is the adiabatic sound speed square. For a > the sound speed square is positive, giving rise to a 
(non-observed) oscillatory behavior. For a < (our case) this square is negative which results in 
instabilities which are not observed either. 

For the viscous fluid the pressure perturbations are qualitatively different. In particular, the 
propagation of perturbations is no longer given by the adiabatic sound speed. The perturbations 
are non-adiabatic. The non-adiabaticity is characterized by 



where 5H = 5 (i%i) /3. In general, the right hand side of (|92p does not vanish (For a similar 
situation see [^g])- ^ * s this non-adiabatic character of the perturbations which makes the power 
spectrum well behaved on small scales in contrast to the GCG case. A more detailed analysis of 
this feature will be given elsewhere. 

VI. CONCLUSIONS 

We presented a two-component cosmological model in which one component represents the dark 
sector, the other one pressureless baryons. The dark sector is described by a single bulk viscous 
fluid with a constant bulk viscosity coefficient £ = £o =constant. A dynamical system analysis was 
used to show that the assumption of a constant £ does not seem to be too restrictive and that this 
model embodies quite general features and does not represent a very particular configuration. In 
the homogeneous and isotropic background the total energy density (i.e., including the baryons) is 
equivalent to that of a GCG with a = —1/2. But while the perturbation dynamics of a GCG is 
entirely adiabatic, the present model exhibits non-adiabatic features. As a consequence, it can avoid 
the problems that usually plague such type of unified dark matter/dark energy models, namely 
the appearance of (non-observed) small scale oscillations or instabilities. The predictions of our 
model were compared with SNe la data and with the 2dFGRS results. For the comparison with 
SNe la data we employed the so-called gold sample of 182 good quality high redshift supernovae. 




(BV model) 



(92) 
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Our Bayesian statistics analysis leads to the following parameter evaluation at 2a level: qo = 
-0.357^;^, n bQ = 0.055+°;°^ , H = 62.34+^ km/s/A^c and t Q = 15.85+°;^ Gy. The main 
differences in comparison with a similar analysis for the ACDM and the GCG models [38j] are: the 
age of the universe is considerably larger; the deceleration parameter is larger (smaller in absolute 
value); the baryon density is only marginally compatible with the primordial nucleosynthesis result. 
A remarkable feature is the (almost) perfect Gaussian distribution for all parameters with a quite 
small dispersion. 

To compare the results of our perturbation analysis with the observed matter power spectrum we 
have used the BBKS transfer function to fix the initial condition at a very high redshift {z ~ 500), 
deep in the matter dominated phase. Since the BBKS transfer function has been conceived for 
a two-component description of the dark sector, dark matter and a cosmological constant, this 
must be seen with some caution. We have argued that, in spite of this problem, the resulting 
power spectrum preserves its shape and can be compared with observations. A good fitting can be 
obtained for values of go which are in agreement with the type la SNe analysis. The best results are 
found for —0.3 > q$ > —0.6. Qualitatively, our analysis confirms previous results for a simplified 
model of the cosmic medium, in which the baryon component was neglected [27]. 
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